Adiabatic quantum pumping in graphene NIS junctions 
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We investigate adiabatic quantum pumping through a normal metal-insulator-superconductor 
(NIS) junction in a monolayer of graphene. The pumped current is generated by periodic modulation 
of two gate voltages, apphed to the insulating and superconducting regions respectively. In the 
bilinear response regime and in the limit of a thin high insulating barrier, we find that the presence of 
the superconductor enhances the pumped current per mode by a factor of 4 at resonance. Compared 
to the pumped current in an analogous semiconductor NIS junction, the resonances have a 7r/2 phase 
difference. We also predict experimentally distinguishable difi'erences between the pumped current 
and the tunneling conductance in graphene NIS junctions. 
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Adiabatic pumping is a transport mechanism in meso- 
and nanoscale devices by which a finite dc current is gen- 
erated in the absence of an applied bias by low-frequency 
periodic modulations of at least two systemparameters 
(typically gate voltages or magnetic fields) [J Q • In or- 
der for electrical transport to be adiabatic, the period 
of the oscillatory driving signals has to be much longer 
than the dwell time Tdweii of the electrons in the system, 
T = 2ttu!~^ ^ Tdweii- Adiabatic quantum pumping Q 
refers to pumping in open systems in which quantum- 
mechanical interference of electron waves occurs. In the 
last decade, many different aspects of quantum pump- 
ing have been investigated in a diverse range of nanode- 
vices, for example charge and spin pumping in quantum 
dots [1], the relation of quantum pumping to geometric 
(Berry) phases [Bj and the role of electron-electron in- 
teractions . Quantum pumped currents have also been 
studied in hybrid systems consisting of normal-metal (N) 
and superconducting (S) parts, such as NS and SNS junc- 
tions [7h11|. Recently, investigations of quantum pump- 



ing in graphene mono- and bilayers have appeared [12l | . 

In this Letter we investigate adiabatic charge pumping 
through a normal metal-insulating-supcrconductor (NIS) 
junction in graphene. The pumped current is generated 
by adiabatic variations of two gate voltages Uo(t) and 
Vo{t) which change, respectively, the Fermi level in the 
superconducting region and the height of the insulating 
tunnel barrier. The central question we aim to answer 
is what the effect of electron-hole (Andreev) reflection 
is on pumped charge currents in graphene. Using the 
scattering matrix formalism, we calculate the adiabati- 
cally pumped current at zero temperature in the linear 
response regime, i.e., for small variations of the pump- 
ing parameters Uq and Vq, and we compare this with the 
pumped current in the absence of the superconducting 
lead. Our main result is that the presence of the super- 
conducting lead enhances the pumped current per mode 
by a factor of 4 and the total pumped current by a fac- 
tor of 3-\/2/2 at the resonant tunneling condition. Off 
resonance, the pumped current is an order of magnitude 




FIG. 1. (a) Sketch of the graphene NIS junction. The variable 
gate voltage Vb creates the insulating barrier I of length d 
and the gate voltage Uo is applied to the superconducting 
electrode (yellow), (b) Schematic of the energy levels in the 
three regions, (c) The area enclosed in the {Uo, Vb) parameter 
space during one pumping cycle. 



smaller than the analogous current in a semiconductor 
NIS junction. We also find that whereas the conductance 
increases with Uq for thin barriers, the pumped current 
decreases with Uq- This difference might be used to dis- 
criminate between conductance and pumped currents in 
graphene NIS junctions. In the last part of the Letter, we 
briefly comment on the pumped current in the so-called 
specular reflection regime (where Aq > Ep with Aq the 
superconductor gap and Ep the Fermi energy) jl3|. 

Consider the geometry depicted in Fig. [1] A ballis- 
tic sheet of graphene in the (x, y)-plane contains a po- 
tential barrier of height Vq and length d and a super- 
conducting contact in the region a; > 0. The barrier 
can be implemented by employing the electric field effect 
14 . isj via a gate voltage and superconductivity can be 
induced in the region a; > via the proximity effect. 
We assume the potential step Vq to be abrupt on both 
sides, which is justified close to the Dirac point where 
the Fermi wavel eng th Af 3> d and which can be realized 
experimentally jl5|. The conductance G{eV) through 
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such a graphene NIS junction has recently been studied, 
both for potential barriers of finite length [l^ and in the 
limit of a thin barrier 17[. Here we investigate the adi- 
abatically pumped current through the latter junction. 
The pumped current is induced by periodic variations of 
Uo{t) = Uo+SUocos{ujt) and Vb(<) = Vo + dVocos{ujt+(f>). 
The total pumped current / into the normal lead (the left 
contact in Fig. [1} can be expressed as an integral over 
the area A that is enclosed in {Uo,Vo) parameter space 
during one period, and is given by the scattering matrix 
expression [8| 

/ = /Ar = |^/" dU„dVo n„^(;7o,Fo) (1) 

-'^SUoSVosiiKj) ^MUo,Vo) (2) 



with 



dUa dVo dUo dVo 



(3) 



Eq. is valid in the bilinear response regime where 
SUq <^ Uo and SVq <C Vq and the integral in Eq. ([T]) be- 
comes independent of the pumping contour. The indices 
a and /3 sum over all modes in the normal lead and S 
denotes the Landaucr-Biittikcr scattering matrix whose 
elements „^ describe the scattering of an electron in 
mode m in lead /3 to a hole in mode n in lead a. 

The low-energy excitations in the NIS junction close to 
the Dirac point K{K') are described by the 4x4 Dirac- 
Bogoliubov-deGennes Hamiltonian [l^ 



H 



na-EF + U{X) 



Ep 



A(x) 

Uix) - Ha 



(4) 



with Ha ~ —iTT-vpiuxdx + sgn(a)aydy), where sgn{a) 
is ± for a = K{K'), vp denotes the Fermi velocity of 
the quasi-particles, and the potential U{x) = —Uq9{x) + 
Vq9{—x)9{x -I- d). The pair potential A(a;) which couples 
the electron and hole excitations has the form Age*'^ and 
we assume that {Ep + Uo) » Aq, the mean- field con- 
dition for superconductivity. Two regimes can be dis- 
tinguished: Ep Aq where the usual retro Andreev 
reflection dominates, and Ep < Aq so ^ Aq, where 



specular Andreev reflection dominates 
sis focuses on the retro reflection regime unless we ex- 
plicitly mention that we study the specular reflection 
regime. The Hamiltonian (j4|) acts on the four-component 
wavefunction -0a = (tpAa, ^Ba, V'Aa' "V'sa)' where A and 
B denote the two nonequivalent sides of the graphene 
unit lattice and a = K'{K) for a = K{K'). Follow- 
ing Ref. [l3|, we also introduce the dimensionless barrier 
strength x = dVo/{1lvp) which allows us to consider the 
limit of a thin barrier where Vq — > oo and d — !> such 
that X remains finite. 




FIG. 2. (a) The pumped current I^'^ (Eq. (O) as a function 
of the potential barrier strength x in a- graphene (g) NIS junc- 
tion, (b) The analogous pumped current I^J''' (Eq. (1101) ') in a 



semiconductor (sc) NIS junction, where \/2mEFd/h — 10. 



After applying continuity of the wavefunctions at the 
boundaries x = —d and a; = 0, one obtains the reflec- 
tion and transmission coefficients of the NIS junction, see 
Eqns. (9) in Ref. [l3|- At the Dirac point e = 0, where 
e denotes the energy of the electrons measured from the 
Fermi level and defining S = Ep/{Ep + Uq), the 
derivatives of the coefficients for normal and Andreev re- 
flection, r and va, with respect to the gate voltages Uq 
and Vq are given by dri^A)/dUn = {dr(^A)/d6){dS/dUo) 
and dr(^A)/dVo = {dr(A)/dx){dxldVa) with 



dr 
dx 
dr* 
'85 
dvA 
dx 

85 



2e*" sin 7 cos a (cos 2x + sin 7 sin a — i cos a sin 2%) 



-le 



(1 -I- sin 7 sin a cos 2x)^ 
' sin a cos a (cos 2% cos a + i sin 2%) 



(1 -I- sin 7 sin a cos 2x) 2 
—2ie~^'^ sin 7 cos 7 sin a cos a sin 2% 
(1 -I- sin7sinQ!Cos2x)2 ' 
— «e"^cosasinQ!(sin7 + sin a cos 2x) 
cos 7(1 -|- sin 7 sin a cos 2x)^ 



(5) 



Here sin 7 = (5 sin a with a the angle of incidence of the 
electron. Substituting Eq. into Eq. © and integrat- 
ing over the angle of incidence yields the pumped current 
at e = 



tNIS 



7r/2 



da- 



mn 7 sm a cos a 
(1 + sin 7 sin a cos 2xY 



(6) 



(7o=0 



4 cos"* X - 12 cos^ X + 8 V2| cos x| - 3 



16\/2| 



cosx|(2cos^ X — 1)'^ 



(7) 



13[. Our analy- where /g_o 



uje-t 



:5Uo5Vo sin0. Eq. ^ is valid for 



Uq = 0. Note from Eq. ^ that there is no contribution to 
I^^^ from normally incident electrons with a = 0. This is 
a display of the Klein tunneling effect where the Andreev 
reflection coefficient |r^(Q; = 0)| = 1 independent of x 
and Uq. 

Figure [2lja) shows the pumped current, Eq. ([7]), as a 
function of x- Notice that Ig'^, just as the conduc- 
tance in this system [l3|, is periodic in x with a pe- 
riod TT. We see that Ig'^ reaches maximum values when 



Vod^ 



l/2)TThvp, where n is an integer. This is the 
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condition for resonant transmission (i.e., r = 0) when the 
conductance of the system reaches a sharply peaked max- 
imum . Due to the sharp changes in the derivative of 
r the pumped current diverges at this point. Expanding 
Eq. ([7|) with respect to x around x 
scales as 



0, we find that I^'^ 



(8) 



where ao = ^^^'§^ ~ 0.014 and ai = 3(1 - ^) 
0.017. The pumped current thus increases with the bar- 
rier height close to x = 0, which is another manifestation 
of the Klein paradox. Next, we analyze the behavior of 
the current with respect to Uq. We expand Eq. ((6)) with 
respect to Uq around C/q = at x = 0, which results in 



tNIS 



/lg,0 ^bo-bi Uo 



(9) 



where bo = ao and bi = « 0.033. Eq. © shows that 
the pumped current decreases with increasing Uq. 

Figure HKb) shows the corresponding pumped current 
I^J^ in a semiconductor NIS junction. We model the 
insulating region again as a barrier of height Vq and 
width d, and define the dimensionless barrier strength 
X = y^2mVod/?i. Solving the Bogoliubov-deGennes equa- 
tion, matching the wavefunction and its derivative at 
the boundaries, obtaining r and and calculating the 
derivatives with respect to the gate voltages Vb and Uo 
yields the pumped current (at the Dirac point e = and 
for C/o = 0); 



NIS 



Isc.O 



tt/2 



da cos a k^d^x^ 



(10) 



(2x(fc2d2 - x^) - {k^d^ + x^) sin2x) 
(2fc2d2x2 cos2 X + (fc^d^ + x'^) sin^ xf ' 



where Isc.o = {uje/Tr){2m'^d'^/K^)SUoSVosm(l) and k = 
{\/2mEF /h) cos a. Eqns. ([7]) and (fTO]) are the main re- 
sults of this Letter. From Eq. (fTO]) we notice that, in 
contrast to graphenc, the normally incident electrons in 
the semiconductor junction do contribute to the pumped 
current, illustrating the absence of Klein tunneling in a 
semiconductor NIS junction. Figure [2Kb) displays ifj^ 
as a function of the barrier strength x- We observe that 
IZ. also oscillates function of x with a period of 
TT. However, the maxima at resonant transmission occur 
when X = nir and are thus shifted by tt/2 with respect to 
the maxima of I^^^. Notice that I^J^ and I^^^ mostly 
have opposite signs and that ifj^ switches sign at several 
points. In addition, the pumped current in a semiconduc- 
tor NIS junction is roughly one order of magnitude larger 
than in graphene [loj . 

An important question arising when thinking about ex- 
perimental detection of pumped currents is how to dis- 
tinguish them from the conductance G in the system. 
In order to answer this question we explicitly compare 
both quantities. The conductance of the NIS junction in 




(b) 



,r/2 



FIG. 3. Ratio of the pumped current (left) and the conduc- 
tance (right) of the NIS junction and the NIN junction in 
graphene versus x for e = and Uo = 0. 



graphene was considered in Rcf. 
e ^ and Uq = 0, 



[17| and is given by, at 



G 



NIS 



Gn 



A{A + 3) + {A^ - 2/1 - 3) arctan {VA) 



where Gn 



4e Ef'W 



A5/2 

(11) 

with w the width of the sample and 
A = cos 2x. Expanding G^^^ for small x around x = 
at [/q = and for small Uo around [/o = at x = 
yields, respectively, 

G^^^/Go^(4-^) + (16-57^)x^ (12) 
iA-Tr) + i2-Tr/2)Uo/EF. (13) 



r^Nis m 

Ug /(-TO 



We find that Gf^^ (4 - 7r)Go for Uq ^ X ^ some- 
what below the ballistic value Go due to mismatch in 
Fermi wavelength in the normal and superconducting 
leads, as mentioned earlier [l3| . Comparing the scaling 
behavior of G^^^ and Ig^^ as a function of x, Eqns. ^ 
and (jl2p . we find that both transport quantities increase 
with increasing x- However when comparing Eqns. (|9]) 
and ([T3| we see that whereas the pumped current I^^^ 



_ decreases with increasing Uq, the conductance G^^^ 



creases with increasing Uq. Intuitively, switching on Uq 
increases the Fermi level mismatch, which increases the 
conductance [20j however decreases the pumped current. 
This difference can be used to discriminate the pumped 
current from the conductance in an actual experiment. 

We now investigate the influence of the superconduct- 
ing lead on the pumped current by comparing I^^^ with 
the pumped current I^^^ through an entirely normal 
NIN junction in graphene. At the Dirac point e = 
and for Uo = the latter is given by 



7-JV/JV 
9 



Ig.O 



(l-|cosxl)^ ^ 
16| cosxl sin^ x 



(14) 



Similarly, the conductance through a NIN junction is 
given by 



G 



NIN 



G, 



sin X — cos^ X arctanh(sin x) 



sin'^ X 

NIS I J NIN . 



(15) 



Figure[3]shows the ratios l^^SjiNiN qNIS jQNiN 
sus X- In both cases, the superconducting lead enhances 
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the transport reaching a maximum at x = 7r/2. For the 
conductance this maximum enhancement is 2 due to the 
contribution of the holes [2lj, while the enhancement of 
the pumped current rises from 2(16 - ll\/2) « 0.89 at 
X = to a factor of 3^2/2 « 2.12 at x = 7r/2. How- 
ever, when comparing I^^^ (Eq. ^) and Ig^'^ per mode 
at X = 7''/2, i.e., before integration over a, we see that 
the superconducting lead enhances the pumped current 
of each mode by a factor of 4. This last result is due 
to both the holes which contribute a factor of 2 and the 
asymmetry of the NIS junction with respect to injection 
of charge carriers which contribute another factor of 2 Q . 

At this point wc briefly mention the behavior of the 
pumped current as a function of an applied bias volt- 
age eV both in the normal (retro) and in the specular 
Andreev reflection regime [2^. For bias voltages below 
the gap < Ao, the pumped current Ig^^ in the retro 
reflection regime decreases from a finite value at eV = 
(see Eq. (O) to zero at eV = Aq. At eV = Aq, the par- 
ticles are fully Andreev reflected (|r^| = 1 independent 
of X and J7o) and therefore Ig^^ = 0. In the specular 
reflection regime [l^ where Ep < Aq and Uq ^ Ep, the 
pumped current exhibits different behavior. First, the 
pumped current, just as the conductance 17|, is insensi- 
tive to X for energies below the gap. The large mismatch 
in Fermi energies of the normal and the superconduct- 
ing leads already acts as a barrier, and as a result the 
addition of another barrier is irrelevant, explaining this 
behavior. Furthermore, the pumped current I^^^ is zero 
for bias voltages equal to the Fermi level eV = Ep, due 
to the absence of Andreev reflection [l3|, and also for 
eV = Aq, see the discussion above. As a flnal remark, the 
pumped current is several orders of magnitudes smaller 
than in the retro reflection regime. 

Finally, we briefly comment on possibilities for experi- 
mental observation of our predictions. Experiments with 
superconducting electrodes on top of graphene in which 
multiple Andreev reflections were observed have already 
been carried out [2^, 24 1. From these experiments we 



can estimate the order of magnitude of the pumped cur- 
rent. Some typical parameters are uj/{2tt) = 5 GHz, 



Ep = 80 meV, 



Vp = 



10 m/s and barrier width 



d = 10 — 20 nm |15l . |25|. For gate voltages on the or- 
der of 10 meV, the pumped current is on the order of 
lOfA far from the resonant tunneling condition, going 
up to 0.1 — 1 pA or higher close to resonance. 

In conclusion, wc have investigated adiabatic quantum 
pumping in a graphene NIS junction, which is generated 
by periodic modulation of the insulating barrier Vq and 
the Fermi level on the superconductor side. We have 
demonstrated that the presence of the superconducting 
lead can enhance the pumped current per mode by a 
factor of 4 (at resonance) and suggested experimentally 
observable differences between the conductance and the 



pumped current in this system. 

This work has been supported by the Netherlands Or- 
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